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Abstract. In this review, recent progress in theoretical models for blazar emission will be summarized. The salient 
features of both leptonic and lepto-hadronic approaches to modeling blazar spectral energy distributions will be 
reviewed. I will present sample modeling results of spectral energy distributions (SEDs) of different types of blazars 
along the blazar sequence, including Fermi high-energy 7-ray data, using both types of models. Special emphasis will 
be placed on the implications of the recent very-high-energy (VHE) 7-ray detections of non-traditional VHE 7-ray 
blazars, including intermediate and low-frequency-peaked BL Lac objects and even fiat-spectrum radio quasars. Due 
to the featureless optical spectra of BL Lac objects, the redshifts of several BL Lacs remain unknown. I will briefly 
discuss possible constraints on their redshift using spectral modeling of their SED including Fermi + ground-based 
VHE 7-ray data. It will be shown that in some cases, spectral modeling with time-independent single-zone models 
alone is not sufficient to constrain models, as both leptonic and lepto-hadronic models are able to provide acceptable 
fits to the overall SED. Subsequently, recent developments of time-dependent and inhomogeneous blazar models will 
be discussed, including detailed numerical simulations as well as a semi-analytical approach to the time-dependent 
radiation signatures of shock-in-jet models. 



1. Introduction 

Blazars (BL Lac objects and 7-ray loud flat spectrum ra- 
dio quasars [FSRQs]) are the most extreme class of ac- 
tive galaxies known. They have been observed at all wave- 
lengths, from radio through VHE 7-rays. The broadband 
continuum SEDs of blazars are dominated by non-thermal 
emission and consist of two distinct, broad components: A 
low-energy component from radio through UV or X-rays, 
and a high-energy component from X-rays to 7-rays (see, 
e.g., Figure [1]). 

Blazars are sub-divided into several types, defined by 
the location of the peak of the low-energy (synchrotron) 
SED component. Low-synchrotron-peaked (LSP) blazars, 
consisting of flat-spectrum radio quasars and low- 
frequency peaked BL Lac objects (LBLs), have their syn- 
chrotron peak in the infrared regime, at Vg £ lO^'* Hz. 
Intermediate-synchrotron-peaked (ISP) blazars, consist- 
ing of LBLs and intermediate BL Lac objects (IBLs) 
have their synchrotron peak at optical - UV frequen- 
cies at 10^* Hz < i/g < 10^^ Hz, while High-synchrotron- 
peaked (HSP) blazars, almost all known to be high- 
frequency-peaked BL Lac objects (HBL), have their syn- 
chrotron peak at X-ray energies with i/g > 10^^ Hz 



(Abdo et al. (2010b)). This sequence had first been iden- 
tified by 'Fossati et al. (1998)' and associated also with 
a trend of overall decreasing bolometric luminosity as 
well as decreasing 7-ray dominance along the sequence 
FSRQ — >■ LBL HBL. According to this classifica- 
tion, the bolometric power output of FSRQs is known 
to be strongly 7-ray dominated, in particular dur- 
ing flaring states, while HBLs are expected to be al- 
ways synchrotron dominated. However, while the overall 



bolometric-luminosity trend still seems to hold, recently, 
even HBLs seem to undergo episodes of strong 7-ray domi- 
nance. A prominent example was the recent joint Fermi + 
HESS observational campaign on the HBL PKS 2155-304 



(Aharonian et al. (2009)) 



Figure [T] shows examples of blazar SEDs along the 
blazar sequence, from the FSRQ 3C279 (a), via the LBL 
BL Lacertae (b) and the IBL 3C 66A (c), to the HBL RGB 
J0710-f591 (d). The sequence of increasing synchrotron 
peak frequency is clearly visible. However, the Fermi spec- 
trum of the LBL BL Lacertae indicates a 7-ray flux clearly 
below the synchrotron level, while the SED of the IBL 
3C 66A is clearly dominated by the Fermi 7-ray flux, in 
contradiction with the traditional blazar sequence. 

The emission from blazars is known to be variable 
at all wavelengths. In particular the high-energy emis- 
sion from blazars can easily vary by more than an 
order of magnitude between different observing epochs 



von Montigny et al. (1995) Mukherjee et al. (1997) 



Mukherjee et al. (1999) ). However, high-energy vari- 
ability has been observed on much shorter time 
scales, in some cases even down to just a few min- 



utes ([Aharonian et al. (2007")| [Albert et al. (2007a) [ ). 



The flux variability of blazars is often accompanied by 
spectral changes. Typically, the variability amplitudes 
are the largest and variability time scales are the shortest 
at the high-frequency ends of the two SED components. 
In HBLs, this refers to the X-ray and VHE 7-ray 
regimes. Such differential spectral variability is sometimes 
associated with inter-band or intra-band time lags as 
well as variability patterns which can be characterized 
as spectral hysteresis in hardness-intensity diagrams 
(e.g. jTakahashi et al. (1996)[ 



Kataoka et al. (2000) 
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Fig. 1. Spectral energy distributions of four sub-classes of blazars: a) the FSRQ 3C279 (from CoUmar et al. (2010) ), 
b) the LBL BL Lacertae, (data from Abdo et al. (2010b)), c) the intermediate BL Lac 3C66A (data from 
Acciari et al. (2010c)|, and d) the HBL RGB J0710+591 (data from [Acciari et al. (2010b)|). In Panel a) (3C279), 



lines are one-zone leptonic model fits to SEDs at various epochs shown in the figure. In all other panels, red lines are 
fits with a leptonic one-zone model; green lines are fits with a one-zone lepto-hadronic model. 



Fossati et al. (2000) [ [Zhang et al. (2002)[ ). However, even 
within the same object this feature tends not to be 
persistent over multiple observations. Also in other types 
of blazars, hints of time lags between different observing 
bands are occasionally found in individual observing cam- 
paigns (e.g., [Bottcher et al. (2007)[ [Koran et al. (2009)[), 



but the search for time-lag patterns persisting throughout 
multiple years has so far remained unsuccessful (see, e.g., 
Hartman et al. (2001) for a systematic search for time 
lags between optical, X-ray and 7-ray emission in the 
quasar 3C279). 

2. Basic features of leptonic and lepto-hadronic 
models 

The high inferred bolometric luminosities, rapid vari- 
ability, and apparent superluminal motions provide com- 
pelling evidence that the nonthermal continuum emis- 
sion of blazars is produced in 1 light day sized emis- 



sion regions, propagating relativistically with velocity /3rc 
along a jet directed at a small angle ^obs with respect 
to our line of sight (for details on the arguments for rel- 
ativistic Doppler boosting, see Schlickeiser (1996) I. Let 



r = (1 — /3p)^^/^ be the bulk Lorentz factor of the emis- 
sion region, then Doppler boosting is determined by the 
Doppler factor D ~ (r[l — /3r cos^obs])^^- Let primes de- 
note quantities in the co-moving frame of the emission 
region, then the observed frequency Vohs is related to the 
emitted frequency through i/obs = D v' /{I + z), where z is 
the redshift of the source, and the energy fluxes are con- 
nected through F""^^^ = F^, . Intrinsic variability on a 
co-moving time scale t'^,^^ will be observed on a time scale 
*var — ^var {'^+z)/D. Using the latter transformation along 
with causality arguments, any observed variability leads 
to an upper limit on the size scale of the emission region 
through R<ct%l^D/[l + z). 

While the electron-synchrotron origin of the low- 
frequency emission is well established, there are two funda- 
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mentally different approaches concerning the high-energy 
emission. If protons are not accelerated to sufficiently 
high energies to reach the threshold for p7 pion pro- 
duction on synchrotron and/or external photons and to 
contribute significantly to high-energy emission through 
proton-synchrotron radiation, the high-energy radiation- 
will be dominated by emission from ultrarelativistic elec- 
trons and/or pairs (leptonic models). In the opposite case, 
the high-energy emission will be dominated by cascades 
initiated by P7 pair and pion production as well as pro- 
ton, TT^, and /i^ synchrotron radiation, while primary 
leptons are still responsible for the low-frequency syn- 
chrotron emission (lepto-hadronic models). The following 
sub-sections provide a brief overview of the main radia- 
tion physics aspects of both leptonic and lepto-hadronic 
models. 



2.1. Leptonic models 

In leptonic models, the high-energy emission is pro- 
duced via Compton upscattering of soft photons 
off the same ultrarelativistic electrons which are 
producing the synchrotron emission. Both the syn- 
chrotron photons produced within the jet (the SSC 
process: [Marscher fc Gear (1985)] [Maraschi et al. (1992)] 
Bloom & Marscher (1996) I, and external photons (the 



EC process) can serve as target photons for Compton scat- 
tering. Possible sources of external seed photons include 
the accretion disk radiation (e.g., Dermer et al. (1992)[ 
Dermer & Schlickeiser (1993)), reprocessed optical - 
UV emission from circumnuclear material (e.g., the 
BLR; [SikoFa et al.(1994)[ iBlandford fc Levinson (1995)| 



Ghisellini fc Madau (1996) 
infrared emission from 



sion 
jet 



Dermer et al. (1997)[ ), 
a dust torus 
synchrotron emis- 
from other (faster/slower) regions of the 
itself ( Georganopoulos fc Kazanas (2003)] 



(Blazejowski et al. (2000) I, or 



Ghiselhni fc Tavecchio (2008) ) 



The relativistic Doppler boosting discussed above al- 
lows one to choose model parameters in a way that the 77 
absorption opacity of the emission region is low through- 
out most of the high-energy spectrum (i.e., low compact- 
ness). However, at the highest photon energies, this effect 
may make a non-negligible contribution to the formation 



of the emerging spectrum ( Aharonian et al. (2008)) and 
re-process some of the radiated power to lower frequen- 
cies. Also the deceleration of the jets may have a signifi- 
cant impact on the observable properties of blazar emis- 
sion through the radiative interaction of emission regions 
with different speed ([Georganopoulos fc Kazan as (2003)[ 
Ghisellini et al. (2005) [) and a varying Doppler factor 
(Bottcher & Principe (2009) ). Varying Doppler factors 



may also be a result of a slight change in the jet orientation 
without a substantial change in speed, e.g., in a helical-jet 



configuration (e.g., Villata & Raiteri (1999) I. In the case 



of ordered magnetic-field structures in the emission re- 
gion, such a helical configuration should have observable 



synchrotron polarization signatures, such as the promi- 
nent polarization-angle swing recently observed in con- 
junction with an optical -I- Fermi 7-ray flare of 3C 279 



(Abdo et al. (2010a) I 



In the most simplistic approaches, the underlying 
lepton (electrons and/or positrons) distribution is ad-hoc 
pre-specified, either as a single or broken power-law 
with a low- and high-energy cut-off. While leptonic 
models under this assumption have been successful in 



modeling blazar SEDs (e.g., Ghisellini et al. (1998) ), 
they lack a self-consistent basis for the shape of the 
electron distribution. A more realistic approach con- 
sists of the self-consistent steady-state solution of the 
Fokker-Planck equation including a physically motivated 
(e.g., from shock-acceleration theory) acceleration of 
particles and all relevant radiative and adiabatic cool- 
ing mechanisms (e.g., Ghisellini & Tavecchio (2009) 
Acciari et al. (2009b)[ [Weidinger et al. (2010)f ^ The 
model described in Acciari et al. (2009b)[ based on the 



time-dependent model of Bottcher fc Chiang (2002) has 



been used to produce the leptonic model fits shown in 

Fig.m 

In order to reproduce not only broadband SEDs, 
but also variability patterns, the time-dependent 
electron dynamics and radiation transfer prob- 
lem has to be solved self-consistently. Such time- 
dependent SSC models have been developed by, e.g.. 



Mastichiadis fc Kirk (1997)| 


Kataoka et al. (2000) 


Li fc Kusunose (2000) 




Sokolov et al. (2004) 


External 



radiation fields have been included in such treatments 
in, e.g., [Sikora et al. (2001)[ [Bottcher fc Chiang (2002) 



Sokolov fc Marscher (2005) 



Leptonic models have generally been very success- 
fully applied to model the SEDs and spectral variabil- 
ity of blazars. The radiative cooling time scales (in the 
observers 's frame) of synchrotron-emitting electrons in a 
typical i3 ~ 1 G magnetic field are of order of several hours 
- 1 d at optical frequencies and <,1 hr in X-rays and 
hence compatible with the observed intra-day variability. 
However, the recent observation of extremely rapid VHE 
7-ray variability on time scales of a few minutes poses se- 
vere problems to simple one-zone leptonic emission mod- 
els. Even with large bulk Lorentz factors of ^ 50, causality 
requires a size of the emitting region that might be smaller 
than the Schwarzschild radius of the central black hole of 
the AGN ("Begelm an et al. (2008)]). A s a possible solution, 
it has been suggested ( Tavecchio fc Ghisellini (2008)) that 
the 7-ray emission region may, in fact, be only a small 
spine of ultrarelativistic plasma within a larger, slower- 
moving jet. Such fast-moving small-scale jets could plau- 
sibly be powered by magnetic reconnection in a Poynting- 



flux dominated jet, as proposed by Giannios et al. (2009) 



2.2. Lepto-hadronic models 

If a significant fraction of the jet power is con- 
verted into the acceleration of relativistic pro- 
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tons in a strongly magnetized environment, 
reaching the threshold for p7 pion production, 
synchrotron-supported pair cascades will develop 
( [Mannheim fc IBiermann (1992)[ [Mannheim "(1993)| ). 



The acceleration of protons to the necessary ultra- 
relativistic energies (£'™^''>,10"'^^ eV) requires high 
magnetic fields of several tens of Gauss to constrain 
the Larmor radius Rl = 3.3 x 10^'"' B^^ i?ig cm, 
where B = 10 Bi G, and Ep = 10^^ Eig eV, to be 
smaller than the size of the emission region, typically 
inferred to be i?;<10^^ cm from the observed variabil- 
ity time scale. In the presence of such high magnetic 
fields, the synchrotron radiation of the primary pro- 
tons 



and 



(Aharonian (2000) Mucke fc Protheroe (2000) I 



of secondary muons 



and 



mesons 



(Rachen & Meszaros (1998) 


Miicke & Protheroe (2000) [ 


Miicke & Protheroe (2001)[ 


Miicke et al. (2003) 1 must be 



taken into account in order to construct a self-consistent 
synchrotron- proton blazar (SPB) model. Electromagnetic 
cascades can be initiated by photons from 7r°-decay ("tt" 



cascade"), electrons from the tt^ 



decay 



("tt^ cascade"), p-synchrotron photons ( "p-synchrotron 
cascade"), and /i-, tt- and ii'-synchrotron photons 
("/x^-synchrotron cascade"). 

Miicke & Protheroe (2001) and Miicke et al. (2003) 
have shown that the "tt" cascades" and "tt^ cascades" 
from ultra-high energy protons generate featureless 7-ray 
spectra, in contrast to "p-synchrotron cascades" and "/i^- 
synchrotron cascades" that produce a two-component 7- 
ray spectrum. In general, direct proton and syn- 
chrotron radiation is mainly responsible for the high en- 
ergy bump in blazars, whereas the low energy bump is 
dominanted by synchrotron radiation from the primary 
e~, with a contribution from secondary electrons. 

Hadronic blazar models have so far been very diffi- 
cult to investigate in a time-dependent way because of 
the very time-consuming nature of the required Montc- 
Carlo cascade simulations. In general, it appears that it 
is difficult to reconcile very rapid high-energy variability 
with the radiative cooling time scales of protons, e.g., due 
to synchrotron emission, which is t°y^ = 4.5 x 10^ (1 -I- 
z) D^^ B:['^ E:[g^ s ( [Aharonian (2000)[ ), i.e., of the order 
of several days for ~ 10 G magnetic fields and typical 
Doppler factors D = 10Z?i. However, rapid variability on 
time scales shorter than the proton cooling time scale may 
be caused by geometrical effects. 

In order to avoid time-consuming Monte-Carlo sim- 
ulations of the hadronic processes and cascades in- 
volved in lepto-hadronic models, one may consider 
a simplified prescription of the hadronic processes. 
Kelner fc Aharonian (2008)[ have produced analytic fit 
functions to Monte-Carlo generated results of hadronic 



interactions using the SOFIA code (Miicke et al. (2000)). 
Those fits describe the spectra of the final decay prod- 
ucts, such as electrons, positrons, neutrinos, and photons 
from tt'^ decay. The use of those functions requires a prior 
knowledge of the target photon field (in the SPB model 



the electron-synchrotron photon field) as well as the pro- 
ton spectrum. Once the first-generation products are eval- 
uated, one still needs to take into account the effect of 
cascading, as the synchrotron emission from most of the 
electrons (and positrons) as well as tt^ decay 7-rays are 
produced at ^ TeV energies, where the emission region is 
highly opaque to 77 pair production. A quasi-analytical 
description of the cascades can be found as follows. 

We start with the injection rates of first-generation 
high-energy 7-rays, iV°, and pairs, Qeil), from the an- 
alytical fit functions of Kelner & Aharonian (2008) The 
cascades are usually well described as linear cascades so 
that the optical depth for 77 absorption, Tjj{e) can be 
pre-calculated from the low-energy radiation field. Under 
these conditions, the spectrum of escaping (observable) 
photons can be calculated as 



1 



'77L 



(1) 



where 7V°™ has contributions from the first-generation 
high-energy photon spectrum and synchrotron emission 
from secondaries, iV™ = + N^^ . A fairly accurate 
evaluation of the synchrotron spectrum is achieved with a 
single-electron emissivity of the form j^, oc v^^^ e~'^^'^° with 
eo = b-f^, where b = B/B^nt and B„it = 4.4 x lO^^ G. 
The electron distribution, Ne{'y) will be the solution to the 
isotropicFokker-Planck equation in equilibrium {d{.)/dt = 
0): 

^ UNeh]) = Qeh) + Nrh) + Nehr^- (2) 

For the high magnetic fields of -B>,10 G, as required for 
hadronically dominated 7-ray emission, electron cooling 
will be dominated by synchrotron losses, i.e. ,7 = —1^07^ 
with = {4/3) cax ub / (meC^). We are only interested 
in high-energy particles emitting synchrotron radiation 
at least at X-ray energies. For those particles, the syn- 
chrotron cooling time scale is expected to be much shorter 
than the escape time scale so that the escape term in Eq. 
[5] may be neglected. N^'' (7) in Eq. [5] is the rate of par- 
ticle injection due to 77 absorption, to be evaluated self- 
consistently with the radiation field. In the 77 absorption 
of a high-energy photon of energy e, one of the produced 
particles will assume the major fraction, fj of the pho- 
ton energy. Hence, an electron/positron pair with energies 
71 = e and 72 = (1 — f-y) e is produced. Furthermore re- 
alizing that every non-escaping photon (according to EqjT]) 
will produce an electron/positron pair, we can write the 
pair production rate as 

Nrh) = /abs(ei) {K+N7,)+Us{e2) {K+K!) (3) 

where ei = l/ff, £2 = 7/(1 - f-f) and /abs(e) = 1 - 
(1 — e~'^T^['^l)/T-y-y[e]. With this approximation, we find an 
implicit solution to Equation [2J 



(4) 
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The solution ^ is implicit in the sense that the par- 
ticle spectrum Neij) occurs on both sides of the equa- 
tion as depends on the synchrotron emissivity, which 
requires knowledge of Ne{"f), where pairs at energies of 
71 = Vl/Uy^ and 72 = \/7/([l - A] b) provide the ma- 
jority of the radiative output relevant for pair produc- 
tion at energy 7. However, for practical applications, one 
may use the fact that generally, 7, the argument on the 
l.h.s., is much smaller than 71,2- This condition is ful- 
filled if there is no pion-induced pair injection at ener- 
gies above 7crit = 4.4 x lO^^ ([1 - /^] Bg)"^ or £;e^crit = 
2.3 X lO^s ([1 - f-y] Bg)~^ eV. In a usual synchrotron- 
proton-blazar model setup, no substantial pair injection 
above -Eccrit is expected. Therefore, Eq. 2] may be evalu- 
ated progressively, starting at the highest pair energies for 
which Qoij) 7^ or iV°^ , 7^ 0, and then using the solu- 
tion for Ne{"f) for large 7as one progresses towards lower 
values of 7. 

Once the equilibrium pair distribution -/Ve(7) is known, 
it can be used to evaluate the synchrotron emissivity and 
hence, using Eq. [1] the observable photon spectrum. 

Example model fits to several blazar SEDs us- 
ing a simplified lepto-hadronic model based on the 



3C279 



Kelner & Aharonian (2008) fit functions and the above 
cascade description arc shown in Fig [U b) - d) and 
compared to leptonic models of the same SEDs. As the 
low-frequency component is electron-synchrotron emission 
from primary electrons, it is not surprising that virtu- 
ally identical fits to the synchrotron component can be 
provided in both types of models. In the high-frequency 
component, strongly peaked spectral shapes, as, e.g., in 
3C 66A and RGB J0710+591 require a strong proton- 
synchrotron dominance with the cascading of higher- 
energy TeV) emission only making a minor contri- 
bution to the high-energy emission. This, in fact, makes 
it difficult to achieve a substantial extension of the escap- 
ing high-energy emission into the > 100 GeV VHE 7-ray 
regime. In objects with a smoother high-energy SED, e.g., 
BL Lacertae in Fig.[T]D, a substantially larger contribution 
from cascade emission (and leptonic SSC emission) is al- 
lowed to account for a relatively high level of hard X-ray 
/ soft 7-ray emission. This also allows for a substantial 
extension of the 7-ray spectrum into the VHE regime. 

3. Spectral fits along the blazar sequence 

In the framework of leptonic models, the blazar sequence 
FSRQ LBL IBL HBL is often modeled through 
a decreasing contribution of external radiation fields to 
radiative cooling of electrons and production of high- 



energy emission (Ghisellini et al. (1998)). In this sense, 
HBLs have traditionally been well represented by pure 
synchrotron-self-Compton models, while FSRQs often re- 
quire a substantial EC component. This interpretation 
is consistent with the observed strong emission lines in 
FSRQs, which are absent in BL Lac objects. At the same 
time, the denser circumnuclear environment in quasars 
might also lead to a higher accretion rate and hence a more 
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Fig. 2. Fits to the SED of 3C279 during the time of 



the MAGIC detection ( [Albert et al. (2008)[ ). The heavy 
solid curve is an attempt of a leptonic model with EC- 
dominated 7-ray emission, while the dotted curve is an 
SSC fit to the X-ray through 7-ray spectrum. The dashed 
and dot-dashed curves are hadronic model fits. 



powerful jet, consistent with the overall trend of bolomet- 
ric luminosities along the blazar sequence. This may even 
be related to an evolutionary sequence from FSRQs to 
HBLs governed by the gradual depletion of the circumnu- 



clear environment (Bottcher & Dermer (2002) I. 

However, in this interpretation, it would be expected 
that mostly HBLs (and maybe IBLs) should be detectable 
as emitters of VHE 7-rays since in LBLs and FSRQs, 
electrons are not expected to reach TeV energies. This 
appears to contradict the recent VHE 7-ray detections 
of lower- frequency peaked objects such as W Comae 



( [Acciari et al. (2008)| ), 3C66A ( Acciari et al. (2009a) ) 
PKS 1424+240 (Acciari et al. (2010a)) 



BL Lacertae ([Albert et al. (2007b)f ; S5 0716-^714 
( [Anderhub et al. (209)[ and ewn the FSRQs 
3C 279 ([Albert et al. (2008)[) and PKS 1510-089 



(Wagner fc Behera (2010)) 



The overall SEDs of IBLs detected by VERITAS could 
still be fit satisfactorily with a purely leptonic model. 
Fitting the SEDs of the IBLs 3C66A and W Comae 
with a pure SSC model, while formally possible, would 
require rather extreme parameters. In particular, mag- 
netic fields several orders of magnitude below equiparti- 
tion would be needed, which might pose a severe prob- 
lem for jet collimation. Much more natural fit param- 
eters can be adopted when including an EC compo- 
nent with an infrared radiation field as target photons 



In 



(Acciari et al. (2009b) Acciari et al. (2010c) ) 



Bottcher et al. (2009) it has been demonstrated 



that the VHE 7-ray detection of the FSRQ 3C 279 poses 
severe problems for any variation of single-zone leptonic 
jet model. The SED with the MAGIC points is shown in 
Fig. 121 along with attempted leptonic as well as hadronic 



6 



M. Bottcher: Models for the Spectral Enerpv nistrihntinns and Variahilitv nf Rlazars 



model fits. The fundamental problem for the leptonic 
model interpretation is the extremely wide separation be- 
tween the synchrotron and 7-ray peak frequencies. In a 
single-zone SSC interpretation this would require a very 
high Lorentz factor of electrons at the peak of the electron 
distribution and, in turn, an extremely low magnetic field. 
In an EC interpretation (the solid curve in Fig. [2]) the SSC 
component, usually dominating the X-ray emission in lep- 
tonic fits to FSRQs like 3C 279, is too far suppressed to 
model the simultaneous RXTE spectrum. As an alterna- 
tive, the optical emission could be produced in a separate 
emission component. A pure SSC fit to the X-ray and 7- 
ray component (the short-dashed curve in Fig. [5]) is tech- 
nically possible, but also requires a far sub-equipartition 
magnetic field. Much more natural parameters could be 
achieved in a fit with the synchrotron-proton-blazar model 



PKS 1424+240 

June 2009 



( Milcke et al. (2003) I , as illustrated by the long-dashed 



and dot-dashed model curves in Fig. O 



4. Redshift constraints from blazar SED modeling 

BL Lac objects are defined by the absence of broad emis- 
sion lines in their optical spectra. This makes the deter- 
mination of their redshifts, and hence distances, very dif- 
ficult, and in many cases even impossible. Photometric 
redshifts, estimated under the assumption that the host 
galaxies of BL Lac objects are nearly standard candles 
(e.g.. 



Sbarufatti et al. (2005)1, are still highly uncertain. 



and optical spectroscopy often yields only lower limits on 
redshifts (e.g. 



Finke et al. (2007)). 



An alternative method of estimating BL Lac red- 
shifts exploits the fact that high-energy and VHE 
7-rays are absorbed through 77 pair production 
on the Extragalactic Background Light (EBL), see. 



e-g-, 


Dwek & Krennrich (2005)| 


Stecker et al. (2006) 


Franceschini et al. (2008) Gilmore et al. (2009) 


Finke et al. (2010) 


Substantial progress has been 



made in the last few years so that the accepted models of 
the EBL spectrum are now converging to an EBL level 
near the lower limit set by direct galaxy counts (e.g.. 



Aharonian et al. (2006) Abdo et al. (2010c) I 



Exploiting this improved knowledge of the EBL spec- 
trum and intensity, limits on VHE 7-ray blazars can now 
be set from the expectation that the intrinsic (unabsorbed 
by the EBL) HE to VHE 7-ray spectra of blazars are not 
becoming harder towards higher energies. A robust upper 
limit to the source redshift can therefore be found by con- 
sidering a straight power-law extrapolation of the best-fit 
Fermi power-law spectrum into the VHE 7-ray regime, 
and attributing the change in spectral slope towards the 
measured VHE spectrum purely to EBL absorption (e.g., 
Georganopoulos et al. (2010) Acciari et al. (2010a)). A 



somewhat more realistic approach has been suggested by 
Prandini et al. (2010) taking into account the expected 



intrinsic softening of the VHE 7-ray spectrum, based on a 
sample of blazars with known redshifts detected both by 
Fermi and ground-based VHE 7-ray observatories. 
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Fig. 3. Fits to the SED of PKS 1424-^240, assuming a 



range of plausible redshifts (from Acciari et al. (2010a)). 
All SEDs have been corrected for 77 absorption by the 



EBL using the model of Gilmore et al. (2009) The inset 



shows a zoom-in on the VHE 7-ray spectrum. The spectral 
fits favour a source redshift of 2; < 0.4. 



If the SED of a blazar with unknown redshift is well 
sampled, the EBL absorption effect can be exploited in 
an even more sophisticated way by constraining models 
through the SED from the optical to GeV 7-rays and 
varying the redshift to yield the best fit to the mea- 
sured VHE 7-ray spectrum. This procedure has been car- 
ried out in detail for the recent VERITAS detections 
of PKS 1424-F240 ([Acciari et al. (2010a)[ ) and 3C66A 
([Acciari et al. (2010c)[). 



In the case of PKS 1424+240, only very unreliable red- 
shift estimates were available. Fig. |3| shows our suite of 
models for a wide range of plausible redshifts. While a 
straight power-law extrapolation of the Fermi spectrum 
yields an upper limit of z < 0.6, our spectral model (lep- 
tonic SSC) clearly favors a redshift of z < 0.4. The red- 
shift of 3C 66A is usually quoted as z — 0.444, but is 



actually highly uncertain as well (Bramel et al. (2005)). 
A scan through redshifts yielded a preferred redshift of 
z ^ 0.2 - 0.3 for this object, in good agreement with 



the estimate found by Prandini et al. (2010) The leptonic 
model fit to 3C66A shown in Fig. [TJ; is based on a redshift 
of z = 0.3. 



5. Inhomogeneous jet models 

The complicated and often inconsistent variability fea- 
tures found in blazars provide a strong motivation 
to investigate jet models beyond a simple, spheri- 
cal, one-zone geometry. The idea behind phenomeno- 
logical multi-zone models like the spine-sheath model 



of Tavecchio & Ghisellini (2008) or the decelerating-jet 
model of Georganopoulos & Kazanas (2003) was that dif- 
ferential relativistic motion between various emission 
zones will lead to Doppler boosting of one zone's emis- 
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Fig. 4. Dependence of the SED characteristics of time- 
averaged spectra from the intermal-shock model, on the 
external radiation energy density Ucxt- 
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Fig. 5. Dependence on inter-band cross correlations and 
time lags on the external radiation energy density. 



sion into the rest frame of another zone. This can reduce 
the requirements of extreme bulk Lorentz (and Doppler) 
factors inferred from simple one-zone letponic modeling of 
rapidly variable VHE 7-ray blazars, such as Mrk 501 or 
PKS 2155-304, and has led to successful model fits to the 
SEDs of those sources with much more reasonable model 
parameters. 

However, the models mentioned above do not treat 
the radiation transport and electron dynamics in a time- 
dependent way and do therefore not make any robust 
predictions concerning variability and inter-band cross- 
correlations and time lags. In order to address those issues, 
much work has recently been devoted to the investigation 
of the radiative and timing signatures of shock-in-jet mod- 
els, which will be summarized in the following sub-section. 

5.1. Shock-in-jet models 

Early versions of shock-in-jet models were developed 
with focus on explaining radio spectra of extragalactic 
jets, e.g., by Marscher fc Gear (1985), Their application 
to high-energy spectra of blazars was proposed by 



Graff et al. (2008) ). As the current detailed internal-shock 



Spada et al. (2001) Detailed treatments of the electron 



energization and dynamics and the radiation transfer in a 
standing shock (Mach disk) in a blazar jet were developed 
by ISokolov et al (2004)"] [Sokolov k Marscher (2005)] 
Graff et al. (2008) The internal-shock model discussed 



Mimica et al. (2004)[ |Joshi (2009)] assumes that the 



central engine is intermittently ejecting shells of relativis- 
tic plasma at varying speeds, which subsequently collide. 
Such models have had remarkable success in explaining 
SEDs and time lag features of generic blazars. 

The realistic treatment of radiation transfer in an 
internal-shock model for a blazar requires the time- 
dependent evaluation of retarded radiation fields orig- 
inating from all parts of the shocked regions of the 
jet. The model system is therefore highly non-linear 
and can generally only be solved using numerical simu- 
lations (e.g., [Sokolov et al. (2004) [ [Mimica et al. (2004) [ 



models employ either full expressions or accurate ap- 
proximations to the full emissivities of synchrotron and 
Compton emission, a complete simulation of the time- 
dependent spectra and light curves istime-consuming and 
does therefore generally not allow to efficiently explore a 
large parameter space. General patterns of the SED, light 
curves and expected time lags between different wave- 
length bands have been demonstrated for very specific, 
but observationally very poorly constrained, sets of pa- 
rameters. 

To remedy this aspect, Bottcher & Dermer (2010) 
have developed a semi-analytical internal-shock model for 
blazars. In this model, the time- and space-dependent elec- 
tron spectra, affected by shock acceleration behind the for- 
ward and reverse shocks, and subsequent radiative cool- 
ing, is calculated fully analytically. Taking into account 
all light travel time effects, the observed synchrotron and 
external-Compton spectra are also evaluated fully analyti- 
cally, using a (^-function approximation to the emissivities. 
The evaluation of the SSC emission still requires a two- 
dimensional numerical integration. 

This semi-analytical model allowed the authors to ef- 
ficiently scan a substantial region of parameter space and 
discuss the dependence on the characteristics of time- 
averaged SEDs, as well as cross-band correlations and time 
lags. As an example. Fig. [H shows the dependence of the 
SED characteristics on the external radiation field energy 
density, Uext- In the classical interpretation of the blazar 
sequence, an increasing Uext corresponds to a transition 
from BL Lac spectral characteristics to FSRQ-like charac- 
teristics. The decreasing synchrotron peak frequency and 
increasing Compton dominance found in the parameter 
study are reproducing this effect. 

From the internal-shock simulations, energy- 
dependent light curves could be extracted. Using 
the standard Discrete Correlation Function (DCF: 



Edelson & Krolik (1988) ) analysis, inter-band time lags 



could be extracted for any set of parameters. Fig. [5] shows 
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the dependence of the inter-band time lags between 
optical, X-ray and Fermi 7-ray light curves (top panel) 
and the quality of the cross correlation, characterized by 
the peak value of the DCF (bottom panel) as a function 
of the external radiation energy density. 

One of the most remarkable results of this study was 
that only slight changes in physical parameters can lead to 
substantial changes of the inter-band time lags and even a 
reversal of the sign of the lags. This may explain the lack 
of consistency of lags even within the same source. 
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